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Abstract 

We calculate the /-electron spectral function using a Keldysh formalism for the Falicov-Kimball model in infinite 
dimensions. We study the region close to the quantum critical point on both the hypercubic and Bethe lattices. 
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The Falicov-Kimball model [1] involves conduction 
electrons, which are free to move through the lattice, 
and /-electrons which are immobile. The two electrons 
interact with each other by a Coulomb interaction (of 
strength U) when they are located at the same lattice 
site. The Hamiltonian is (at half filling) 

H = - ^ ^ (cjc, +c]c.) + uY, cUfU^ 



u 



(1) 



where c] (c;) creates (destroys) a conduction electron 
at site i, f} (fi) creates (destroys) a localized electron 
at site i, and t* is the hopping integral [2]. The sym- 
bol Z represents the number of nearest neighbors, and 
(ij) denotes a sum over all nearest neighbor pairs. The 
formalism for solving the conduction-electron Green's 
function was worked out by Brandt and Mielsch [3] . 

Brandt and Urbanek [4] describe how to calculate the 
/-electron spectral function using the Keldysh tech- 
nique to directly determine the greater Green's func- 
tion (along the Keldysh contour) and then Fourier 
transforming to real frequency (see also [5] for a review 
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Fig. 1. Keldysh contour for evaluating the /-electron Green's 
function at time t. The contour runs from to t, then back 
from t to and finally goes along the imaginary axis down to 
— i/3. When we discrctize the matrix operator over the Keldysh 
contour, we evaluate the integrals via a rectangular (midpoint) 
summation. We typically use no more than 2500 time steps on 
the contour. 

of the formalism). The greater Green's function is de- 
fined to be 

G>{t) = -Tr{e-'^^--S.(A)/(t)/t(0))/2,™p, (2) 

with f{t) = exp{itTi.imp)f exp{~itTli,np) and the evo- 
lution operator is given by 



Sc{\) = Tcexp 



dt / dt'c\t)\c(t,t')c{t') 



(3) 



The subscript imp denotes the use of the impurity 
Hamiltonian (no hopping) with the evolution operator 
corresponding to the dynamical mean field [3] X{io) (the 
dynamical mean field mimics the hopping of conduc- 
tion electrons onto and off of a given site; it is adjusted 
so that the impurity conduction-electron Green's func- 
tion equals the local lattice Green's function). The 
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time-ordering is along the Keldysh contour (see Fig. 1), 
and the contour-ordered dynamical mean field is found 
from a Fourier transform of A(a') 

oo 

Ac(f, ?) = — — / du!lmX{ui) exp[—iLj(t — ¥)] 
V J 

~ oc 

x[fpD{io)-e,{t^i% (4) 

where fpoiu)) = + exp(/3a;)] is the Fermi-Dirac 
distribution and 9c{t — F) = if F is in front of t on 
the contour c and 1 if it is behind. The Green's func- 
tion can be solved by directly evaluating the Feynman 
path integral over the contour, which reduces to cal- 
culating the determinant of a continuous matrix oper- 
ator. This operator is discretized, and a conventional 
determinant is evaluated when wc compute results nu- 
merically. Once (t) is found, then the DOS (at half 
filling) is determined by 



oo 

Af{u;) = -^ j dtRe{G^(t)}cos(wt). 



(5) 



We work with two different infinite coordination 
number lattices here: (i) the Bethe lattice, which has 
a noninteracting DOS that is a semicircle pb{ijj) = 
\/At*'^ — up' I'livt*'^ and (ii) the hypercubic lattice, 
which has a noninteracting DOS that is a Gaussian 
Ph{oj) = &q){—ijj^/t*^)/^t*. All energies are mea- 
sured in units of t* . 
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Fig. 2. F-electron DOS for the Bethe lattice. The different 
thickness lines correspond to different temperatures. The con- 
duction electron DOS is the dashed line. 

In Fig. 2, wo plot the DOS at the critical value of U 
for the Mott transition on the Bethe lattice (U = 2), 
and in Fig. 3 just on the insulating side of the transi- 
tion on the hypercubic lattice (U = 1.5). Also included 
is a plot of the conduction-electron DOS, which is in- 
dependent of temperature [6] . Note how the /-electron 
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Fig. 3. Similar plot of the /-electron DOS for the hypercubic 

lattice. 

spectral function has an interesting temperature evo- 
lution, where it develops a pseudogap as T is reduced. 
Calculations become more difficult at lower temper- 
atures [7] (because of discretization errors along the 
Keldysh contour), and rapidly exhaust our computer 
tional resources. As a check on the accuracy, we use 
the DOS to calculate the Matsubara frequency Green's 
functions and compare them to results calculated di- 
rectly on the imaginary axis, and wc compare the first 
three moments to the exact results for those moments. 
In general, we need to extrapolate the discretization 
size to zero to get accurate results, and most reported 
DOS here are accurate to at least 1% in all of the mo- 
ments and 0.1% for all Matsubara frequency Green's 
functions. A long paper evaluating different values of 
U is in preparation [7] . 
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